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Due to the rapid longitudinal expansion of the quark-gluon plasma created in relativistic heavy 
ion collisions, potentially large local rest frame momentum-space anisotropics are generated. The 
magnitude of these momentum-space anisotropics can be so large as to violate the central as- 
sumption of canonical viscous hydrodynamical treatments which linearize around an isotropic 
background. In order to better describe the early-time dynamics of the quark gluon plasma, one 
can consider instead expanding around a locally anisotropic background which results in a dynam- 
ical framework called anisotropic hydrodynamics. In this proceedings contribution we review the 
basic concepts of the anisotropic hydrodynamics framework presenting viewpoints from both the 
phenomenological and microscopic points of view. 
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1. Introduction 



The space-time evolution of matter created in relativistic heavy-ion collisions seems to be very 
well described by relativistic viscous hydrodynamics [1-13]. Nevertheless, in the early stages of the 
collisions, due to large flow gradients the viscous corrections to the perfect-fluid energy-momentum 
tensor become very large and the system becomes highly anisotropic in momentum space (in the 
local rest frame). Typically, the transverse pressure, P±, is much larger than the longitudinal pres- 
sure, P||, where the directions are specified with respect to the beam axis. Generally speaking, large 
momentum-space anisotropies pose a problem for 2nd-order viscous hydrodynamics, since it relies 
on a linearization around an isotropic background. It has been shown that large linear corrections 
generate unphysical results such as negative particle pressures, negative one-particle distribution 
functions, etc. [14]. 

A similar aspect of the same problem is that microscopic models of the early stages of rel- 
ativistic heavy-ion collisions suggest that the produced system is highly anisotropic [15]. In the 
theory of the Color Glass Condensate (CGC), at very early proper-times, t -C 1/2^, where Qs is 
the saturation scale, the classical gluon fields lead to an energy-momentum tensor which has the 
form T^'■'^ = diag(e, e, e, — e) [16, 17]. This form implies a negative value of the longitudinal pres- 
sure with the transverse pressure equal to the energy density. At later proper times, T ^ l/2.v» 
both the analytical perturbative approaches [18] and the full numerical simulations [19] lead to the 
form T^^^ = diag(£,e/2,£/2,0). This form implies that the longitudinal pressure is zero. Conse- 
quently, proper matching of the results of the microscopic models with the hydrodynamic descrip- 
tion (where the energy-momentum tensor should be close to the isotropic form) is not straightfor- 
ward. 

In addition to the issues faced at early times, one finds that as one approaches the transverse 
edge of the interaction region, corrections to ideal hydrodynamics can be large at all times. As 
a result, one expects large momentum-space anisotropies to be present at the edges [14] and, as 
a result, one is motivated to find better approximation schemes to describe the dynamics in these 
regions. Finally, and perhaps most importantly, physically one expects entropy production to vanish 
in two limits: the ideal hydrodynamical limit (vanishing shear viscosity) and the free streaming 
limit (infinite shear viscosity); however, within viscous hydrodynamics, entropy production is a 
monotonically increasing function of the shear viscosity. In the large shear viscosity limit, as 
one might expect, viscous hydrodynamics becomes a poor approximation and one is, once again, 
motivated to find an alternative framework. 

The difficulties eluded to above motivated the development of a reorganization of viscous 
hydrodynamics in which one incorporates the possibility of large momentum-space anisotropies at 
leading order. The inclusion of large anisotropies also allows for direct matching with the theories 
such as CGC. This framework has been dubbed anisotropic hydrodynamics [20-26]. Below we 
will briefly review the basic ideas underpinning the framework and direct the reader to various 
references for more details. A more self-contained presentation of the concepts of anisotropic 
hydrodynamics can be found, e.g., in Ref. [27]. 
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2. Anisotropic hydrodynamics 

The framework of anisotropic hydrodynamics may be introduced in (at least) two ways. The 
first method is phenomenological and refers directly to tensor structures describing an anisotropic 
fluid [20, 22, 24, 26]. The second method is microscopic; one uses the Boltzman equation for the 
parton (gluon) distribution function and makes an expansion around the anisotropic background 
[21,23,25]. In this Section we outline the first method. The microscopic approach will be dis- 
cussed in the next Section. 

The dynamic equations of anisotropic hydrodynamics follow from the equations 

d^T"" = 0, (2.1) 
^^^a^ = I, (2.2) 

which express the energy-momentum conservation and entropy production laws, respectively. The 
energy-momentum tensor T^^^ has the form 



(2.3) 



We note that Eq. ( ^j| ) reproduces the energy-momentum tensor of the perfect fluid if the two pres- 
sures are equal. Similarly, in the case £ = the entropy production law in {12) is reduced to 
the entropy conservation law. The four- vector describes the flow of matter and defines 
the beam (z) axis. In the general (3-1-1 )-dimensional case, we use the following parameterizations: 
= (MoCoshT>,M.v,Mv,Mosinht>) and V = (sinh t?, 0,0, cosh tJ), where Ux and Uy are the trans- 
verse components of the four-velocity field and ■& is the longitudinal fluid rapidity. The entropy 
flux equals oU^ , where a is the non-equilibrium entropy density. 

One can show [20] that instead of P\\ and Pj^ it is more convenient to use the entropy density 
a and the anisotropy parameter x as two independent variables (one may use the approximation 
P\\/Pi_ ~ x^^/'^). Similarly to standard hydrodynamics with vanishing baryon chemical potential, 
the energy density e, the entropy density a, and the anisotropy parameter x are related through the 
generalized equation of state [22] 



e{x,o) = Cqgp (a)r(x), 
P^{x,o) = Pqgp(cy) [r(x)+3xr'(x) 
P\\{x,o) = Pqgp(a) [r{x)-6xr'{x) 



(2.4) 



where the functions £qgp and Pqgp correspond to a realistic isotropic QGP equation of state, see 
e.g. Ref. [28]. The function r{x) characterizes properties of the fluid which exhibits the pressure 
anisotropy specified by the variable x [20, 29, 30] 



r{x) 



X 3 



1 + 



;carctan \Jx — 1 



(2.5) 



In the isotropic case x= 1, r(l) = 1, / (I) = 0, and Eq. (2.4) is reduced to the standard equation of 
state used in [28]. 
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The function £ appearing in Eq. (2.2) defines tlie entiopy production. In the phenomenological 



approach [20] the following form of E has been used 

V ^ '^eq 

where the time-scale parameter Tgq controls the rate of equilibration. In Refs. [21, 23, 25] a medium 
dependent Teq was used, which was inversely proportional to the typical transverse momentum scale 
in the system. In this case, as we will discuss in some detail below, the late time dynamics is that 
of second order viscous hydrodynamics. If a constant value of Teq is used, the system instead 
approaches perfect fluid behavior for T » Tgq. 

If the generalized equation of state and the entropy production term are defined, the dynam- 
ical equations of anisotropic hydrodynamics become 5 equations for 5 unknown functions: three 
independent components of the fluid four- velocity i/'^, the entropy density a, and the anisotropy 
parameter x. In each case one has to additionally specify the value of the relaxation time or its 
dependence on other thermodynamics-like variables such as a or x. 



3. Microscopic approach 



In the microscopic approach [21, 23, 25] one considers locally anisotropic momentum distribu- 
tions of partons that are characterized by two different scales: X± and They may be interpreted 
as the transverse and longitudinal temperatures {LRF stands below for the local rest frame) 



flRF — f 



Pi- \P\\ 



(3.1) 



In the studies of anisotropic quark-gluon plasma one often uses the Romatschke-Strickland form 
[31] 



/lrf — f 



\ 




f 



1 



p]_+xpj 



(3.2) 



where / is an arbitrary isotropic distribution function which is set by the underlying equilibrium 
distribution function for the particles being considered. The covariant version of (3.2) reads 



/ = / 



1 



{p.ur + ^ip-vy 



(3.3) 



where ^ = l—x. Using Eq. (33) in the standard definitions of the energy-momentum tensor and 
the entropy flux one obtains the expressions for T^^ and a'' of the form introduced in the previous 
Section. 

In order to obtain the dynamical equations in the microscopic approach, one uses the Boltz- 
mann equation. The zeroth moment of the Boltzmann equation, p^d^f = C, with the collision term 
treated in the relaxation-time approximation, C = —p ■ UT {f — fsq), gives 



(3.4) 
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Here F is the inverse relaxation time (F is inversely proportional to Tgq introduced earlier in the 
entropy source, however, the proportionality constant might be different from unity), /gq is the 
background equilibrium distribution function, and n is the parton density (ngq is the corresponding 
equilibrium parton density). We note that for the Romatschke-Strickland form one finds a = 4n. 
The first moment of the Boltzmann equation yields the energy-momentum conservation law 

d^jdPp'p^'f = d^r^=0. (3.5) 

These four equations are satisfied only if the collision term fulfills the condition 

j dPp''C=- j dPp-UpTif- /eq) = 0. (3.6) 

The last equation is the Landau matching condition that specifies the temperature, T , appearing in 
the background equilibrium distributions /eq. If the Landau matching is used to eliminate T, we are 



left with Eqs. ( [3.4| ) and ( [3.5| ) which are 5 equations for 5 unknowns, similarly to the case described 
in the previous Section. 

The structure of the equations of the anisotropic hydrodynamics obtained from the microscopic 
considerations is the same as that used in the phenomenological approach. On the other hand, the 
advantage of the microscopic approach is that the explicit form of the distribution function defines 
the entropy source term £ and the generalized equation of state. If the Romatschke-Strickland form 
is used, the generalized equation of state has the form {1A)-(1S), however, the functions £qgp and 
Pqgp describe a massless parton gas (with Bose or Boltzmann statistics). The phenomenological 
and microscopic approaches have been compared in more detail in [34]. 

Within the microscopic approach one can also perform an explicit matching to Israel-Stewart 
second order viscous hydrodynamics in the limit of small anisotropy. In Ref. [21] it was shown 
that by linearizing the (O-i-l)-dimensional anisotropic hydrodynamics dynamical equations and re- 
quiring that they reproduce second-order viscous hydrodynamics, one can fix the relaxation rate F 



appearing in Eq. (3.4) to be 



^=T' ^^^ = 4^' ^^-^^ 

^Tt ^ ' eq 

which for an ideal equation of state results in 

IT 

r= — , (3.8) 

where f] = T]/5'eq with r] being the shear viscosity, ^eq is the equilibrium entropy density, and T is 
the effective temperature determined via Landau matching in the first moment as described above. 
Two points can be made using this result: (1) the relaxation rate F is proportional to the effective 
local temperature (fourth root of the energy density) and hence in regions of lower temperature, 
the relaxation to isotropy will be slower and (2) the relaxation rate F is inversely proportional to 
the shear viscosity to entropy ratio and hence when the shear viscosity increases, one expects the 
relaxation to isotropy will be slower. Both of these effects are important to understanding the 
dynamics of isotropization in the quark gluon plasma. 
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We want to emphasize that the microscopic approach presented here uses the zeroth and first 
moments of the Boltzmann equation. The zeroth moment in the standard approach to relativistic 
fluid dynamics is used to introduce the particle (or baryon) number conservation law. This is done 
by introducing a chemical potential together with the appropriate Landau matching condition 
which can be used to establish the value of n (in a similar manner as the temperature was defined 
above). In our opinion this procedure is not appropriate for the systems produced in the early 
stages of relativistic heavy-ion collisions. Such systems consist primarily of gluons whose number 
is not conserved. In fact, the radiation of soft gluons leads to the increase of their numbers and 
to entropy production. Hence, the production of gluons is the main physical process which leads 
to equilibration of the system. Interestingly, such a physical picture is qualitatively similar to the 
bottom-up thermalization scenario [35]. 

4. Summary 

In this brief proceedings contribution we have attempted to review the motivation for and the 
foundational ideas behind anisotropic hydrodynamics. The development of this framework has oc- 
curred using two separate methods, dubbed the phenomenological approach and the microscopic 
approach herein. We have reviewed both methods pointing out the fundamental connections that 
exist between the two approaches. In both cases presented, the framework of anisotropic hydrody- 
namics is constructed in such a way that the pressure (momentum) anisotropics of the system are 
built directly into the formalism. In our opinion this framework forms an alternative to the standard 
viscous-hydrodynamics approach where one starts with the isotropic distribution and introduces 
potentially large corrections in the expansion. 

The use of anisotropic hydrodynamics may help us to circumvent the early thermalization 
puzzle. It turns out that all soft hadronic observables may be successfully reproduced in the ap- 
proach where the system is initially highly anisotropic [26]. It turns out that the effects of the initial 
anisotropy may be compensated by a change of the initial conditions. This points to a possible uni- 
versality of the results for flow with respect to initial condition variation (at least with respect to 
the initial momentum-space anisotropy which is assumed). 

There have been many other recent developments which we were not able to go into much 
detail about herein. The approach presented has been generalized to incorporate the presence of 
dynamical longitudinal chromoelectric oscillations [32] and to include a mixture of two interacting 
anisotropic (quark and gluon) fluids which are characterized by different anisotropy parameters and 
equilibration times [33]. The work on mixtures of coupled quark and gluon fluids is potentially im- 
portant because it demonstrates how to enforce number conservation in the baryonic (quark) sector 
while allowing number non-conservation in the gauge (gluon) sector. An interesting project for the 
future would be to find more connections between anisotropic hydrodynamics and the studies of 
the anisotropic systems based on the AdS/CFT correspondance [36^1]. 
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